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Lighter-than  air  vehicles  are  an  attractive  solution  for  many  applications  requiring  a 
sustained  airborne  presence.  The  buoyancy  force  provides  an  energy-free  form  of  lift, 
offering  a  non-traditional  approach  to  long-duration  missions  for  which  traditional  aircraft 
are  not  well-suited.  Potential  applications  include  roving  or  hovering  surveillance  and 
communication  utilities  for  both  military  and  commercial  use,  and  a  variety  of  remote¬ 
sensing  instruments  for  the  scientific  community.  In  particular,  the  Missile  Defense  Agency 
plans  to  utilize  unmanned  airships  at  high-altitudes  to  provide  a  long-duration  missile 
defense  presence  around  the  coast-line  of  the  United  States.  Operated  at  70  kft,  each  of 
these  “high  altitude  airships”  will  fly  above  all  regulated  air-traffic  for  several  months  to 
years,  will  reside  in  a  steady  atmospheric  regime,  and  will  utilize  solar  energy  to  provide 
all  required  power. 

Two  key  objectives  for  this  type  of  mission  are  that  the  unmanned  airship  have  excep¬ 
tionally  long  endurance,  and  that  it  operate  with  a  sufficiently  high-level  of  autonomy.  In 
order  to  achieve  these  objectives,  a  robust  guidance  and  control  system  is  required,  capable 
of  auto-piloting  and  controlling  the  airship  under  an  extremely  wide  range  of  atmospheric 
and  wind  conditions.  The  successful  design  of  such  a  system  first  requires  an  accurate 
model  of  airship  dynamics  across  its  expansive  flight  envelope,  and  a  representative  model 
of  the  expected  disturbances. 

The  dynamics  of  an  airship  are  markedly  different  from  traditional  aircraft,  with  signif¬ 
icant  effects  from  added  mass  and  inertia,  and  a  much  higher  sensitivity  to  wind.  In  this 
paper,  a  typical  airship  configuration  is  first  sized  to  meet  energy  balance  and  mass  con¬ 
straints.  The  geometry  of  this  configuration  is  then  used  to  develop  a  general  aerodynamic 
model  for  the  airship.  The  equations  of  motion  with  added  mass  and  inertia  are  developed, 
and  the  open-loop  dynamics  are  analyzed  across  a  range  of  flight  conditions.  Finally,  con¬ 
trol  laws  are  designed  for  a  single  operating  condition,  and  the  closed-loop  performance  is 
presented  across  a  range  of  velocities. 

I.  Introduction 

Lighter-than-air  (LTA)  vehicles  represent  a  unique  and  promising  platform  for  many  applications  that 
involve  a  long-duration  airborne  presence.  Achieving  their  lift  through  buoyancy,  these  vehicles  require 
much  less  power  than  traditional  aircraft.  This  opens  the  door  to  renewable  sources  of  energy,  which  would 
enable  controlled  flight  to  be  maintained  indefinitely.  Potential  applications  include  roving  or  hovering 
surveillance  and  communication  utilities  for  both  military  and  commercial  use,  plus  a  variety  of  remote- 
sensing  instruments  for  the  scientific  community.  In  particular,  the  Missile  Defense  Agency  (MDA)  plans  to 
utilize  unmanned  airships  at  high  altitudes  to  provide  a  long-duration  missile  defense  presence  around  the 
coast-line  of  the  United  States.  Operated  at  70  thousand  feet,  each  of  these  “high  altitude  airships”  (HAA) 
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will  fly  above  all  regulated  air-traffic  for  several  months  to  years,  will  reside  in  a  steady  atmospheric  regime, 
and  will  utilize  solar  energy  to  provide  all  required  power. 

The  HAA  concept  has  several  attractive  features,  making  it  a  promising  solution  for  both  government 
programs  and  commercial  ventures.  The  large  surface  of  the  HAA  lends  itself  to  the  integration  of  solar  cells. 
The  operational  altitude  is  above  all  cloud  formations,  enabling  direct  sunlight  to  reach  the  solar  cells,  and  it 
resides  in  a  calm  portion  of  the  atmosphere  where  the  wind  speeds  are  relatively  low.  The  HAA  concept  has 
several  advantages  over  satellites  as  well.  It  can  maintain  a  geostationary  location  much  closer  to  the  earth, 
matching  both  the  longitude  and  latitude  at  any  point  on  the  globe.  Geostationary  satellites  are  about  1,700 
times  further  away,  and  they  have  limited  coverage  for  extreme  north/south  latitudes.  Continuous  coverage 
could  be  provided  from  low-earth  orbit,  but  would  require  a  constellation  of  satellites.  In  addition,  the  HAA 
has  the  flexibility  to  change  its  target  location  over  time,  and  it  may  be  grounded,  repaired,  and  returned 
to  flight  at  any  time  if  necessary. 

Key  ingredients  for  success  of  the  HAA  concept,  whatever  the  application,  are  that  it  have  exceptionally 
long  endurance,  and  that  it  operate  with  a  sufficiently  high-level  of  autonomy.  In  order  to  achieve  these 
objectives,  a  robust  guidance  and  control  system  is  required,  capable  of  auto-piloting  and  controlling  the 
airship  under  an  extremely  wide  range  of  atmospheric  and  wind  conditions.  The  successful  design  of  such 
a  system  first  requires  an  accurate  model  of  airship  dynamics  across  its  expansive  flight  envelope,  and  a 
representative  model  of  the  expected  disturbances. 

The  dynamics  of  an  airship  are  markedly  different  from  traditional  aircraft,  with  significant  effects  from 
added  mass  and  inertia,  and  a  much  higher  sensitivity  to  wind.  In  this  paper,  a  typical  airship  configuration 
is  first  sized  to  meet  energy  balance  and  mass  constraints.  The  geometry  of  this  configuration  is  then  used 
to  develop  a  general  aerodynamic  model  for  the  airship.  The  equations  of  motion  with  added  mass  and 
inertia  are  developed,  and  the  open-loop  dynamics  are  analyzed  across  a  range  of  flight  conditions.  Finally, 
control  laws  are  designed  for  a  single  operating  condition,  and  the  closed-loop  performance  and  robustness 
is  presented  with  respect  to  a  range  of  velocities. 

A.  Mission  Scenario 

The  mission  scenario  used  in  this  paper  is  summarized  in  Table  1.  It  is  consistent  with  the  mission  scenario 
identified  by  the  MDA. 


Table  1.  Mission  Scenario 


English 

Metric 

Altitude 

70,000  kft 

21.3  km 

Payload 

4,000  lb 

1,814  kg 

Power  to  Payload 

- 

15  kW 

Nominal  Airspeed 

35  knots 

18  m/s 

Maximum  Airspeed 

90  knots 

46  m/s 

II.  High-Level  Airship  Design 

In  general,  airships  fall  into  one  of  two  categories:  pressurized  or  rigid.  With  pressurized  models,  the 
shape  of  the  hull  is  maintained  by  a  small  pressure  differential  between  the  internal  lifting  gas  and  the 
external  atmosphere.  With  rigid  airships,  the  hull  is  formed  by  a  lightweight  skeletal  structure  with  a  fabric 
skin,  while  the  lifting  gas  is  stored  in  internal  “gasbags”.  The  Graf  Zeppelin  II  and  the  Hindenburg  were 
both  rigid  airships.  Pressure  airships  are  sometimes  made  “semi-rigid”  by  adding  a  structural  keel  to  share 
the  bending  loads  generated  from  aerodynamic  forces.  In  this  paper,  a  pressurized  design  is  assumed. 

The  classic  airship  design  consists  of  an  axisynnnetric,  teardrop-shaped  hull  with  a  hanging  empennage  (or 
gondola)  and  tail  fins  for  stability.  For  the  purpose  of  this  analysis,  the  hull  is  modeled  as  two,  axisynnnetric 
ellipsoid  halves.  Each  ellipsoid  has  the  same  semi-minor  axis,  b ,  but  different  semi- major  axes,  aq  and  aq. 
Figure  1  on  the  next  page  shows  a  2-D  view  of  the  double  ellipsoid  model.  More  complex  geometries  have 
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Figure  1.  Double-Ellipsoid  Geometry  for  Airship  Hull 


been  designed  with  the  aid  of  CFD  tools,  with  the  primary  goals  of  reducing  turbulence  and  drag.1  Prior  to 
the  advent  of  CFD,  however,  ellipsoid  geometries  had  historically  been  used  to  study  the  aerodynamics  of 
airships.  The  ellipsoid  configuration  is  therefore  used  in  this  analysis  as  it  captures  the  essential  aerodynamic 
characteristics  of  typical  airships,  while  facilitating  a  mathematical  development  of  the  model. 

The  first  step  in  the  airship  design  process  is  to  size  the  hull.  The  volume  of  the  hull  defines  the  buoyant 
lift  capability  of  the  airship,  and  determines  the  maximum  attainable  altitude.  It  is  first  necessary,  however, 
to  develop  some  relationships  between  mass,  volume,  and  gas  densities,  as  governed  by  aerostatic  principles. 
Next,  with  the  size  and  shape  of  the  airship  defined,  we  may  compute  the  drag  force  experienced  at  a 
particular  airspeed.  This  in  turn  drives  the  thrust  and  power  requirements.  The  required  coverage  area  for 
solar  cells  may  then  be  analyzed,  along  with  the  mass  available  for  structure  and  sub-systems.  The  objective 
of  the  analysis  is  to  ensure  that  the  proposed  configuration  is  feasible. 

A.  Aerostatics  of  Airships 

The  upward  buoyancy  force  generated  by  an  airship  is  equal  to  the  weight  of  the  displaced  air.  This  force  is 
typically  referred  to  as  the  “gross  lift”,  and  is  defined  as: 

Lg  =  VnPa9  (1) 

where  Vjv  is  the  net  volume  of  displaced  air,  pA  is  the  density  of  air,  and  g  is  the  acceleration  due  to  gravity. 
Subtracting  the  weight  of  the  lifting  gas  (Helium),  we  obtain  the  “net  lift”,  L at.  Noting  that  the  volume  of 
displaced  air  is  equal  to  the  volume  occupied  by  the  Helium,  the  equation  for  the  net  lift  is: 

Ln  =  VN  (pA  -  pH)  g  (2) 

This  represents  the  amount  of  lift  available  to  counteract  the  weight  of  the  airship  structure  and  payload.2 

The  density  of  both  the  air  and  Helium  vary  with  altitude.  Assuming  that  both  gases  have  the  same  pres¬ 
sure  and  temperature,  their  densities  change  uniformly  with  altitude.  Although  a  slight  pressure  differential 
is  required,  it  is  small  enough  to  neglect  for  the  purposes  of  this  analysis. 

As  the  airship  rises,  the  density  of  Helium  decreases  along  with  the  atmospheric  density.  Therefore,  since 
the  mass  of  the  Helium  remains  fixed,  the  volume  Vjv  must  increase.  This  variation  in  internal  lifting  gas 
volume  is  achieved  through  ballonets  bags  of  air  inside  the  hull  which  expand  and  contract  to  regulate  the 
internal  pressure  (and  thereby  the  volume).  At  the  launching  altitude  (assume  sea  level),  the  density  is  at 
its  highest  value.  The  ballonets  are  expanded  to  their  maximum  volume,  and  Vn  is  at  a  minimum.  As  the 
airship  begins  to  rise,  the  ambient  density  and  pressure  both  fall,  and  air  is  automatically  ejected  from  the 
ballonets  to  match  the  falling  pressure.  Clearly,  at  some  point  during  the  ascent,  the  ballonets  will  become 
completely  empty.  At  this  point,  no  further  expansion  of  the  lifting  gas  volume  is  possible.  The  net  volume 
has  reached  a  maximum  value,  Vmax .  This  point  is  termed  the  “pressure  altitude” . 

Continued  ascent  causes  a  reduction  in  net  left  as  the  densities  fall  but  the  volume  remains  constant. 
The  pressure  differential  also  increases,  creating  a  superpressure  condition,  which  can  result  in  rupture  of 
the  exterior  skin  if  the  differential  becomes  too  great.  Helium  may  be  vented  to  avoid  rupture,  but  this  is 
extremely  undesirable  in  long-endurance  applications  as  it  reduces  the  available  lift  and  shortens  the  mission 
life.  It  is  therefore  an  important  criteria  in  autonomous  operations  to  keep  the  airship  below  the  pressure 
altitude. 
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It  can  be  shown  that  the  net  lift  is  constant  over  all  altitudes,  up  to  the  pressure  altitude.  This  is  based 
upon  the  assumption  that  the  density  of  the  lifting  gas  changes  at  the  same  rate  as  the  atmospheric  density. 
Let  pao  be  the  density  of  air  at  sea  level,  and  denote  the  density  ratio  as: 

°  =  —  (3) 

PAO 

where  pA  is  the  density  of  air  at  a  given  altitude.  The  density  of  Helium  at  that  altitude  is: 

Ph  =  °Pa  (4) 

Noting  that  the  net  volume  is  equivalent  to  the  mass  of  Helium  divided  by  the  density  of  Helium,  we  may 
rewrite  Eq.  2  on  the  preceding  page  as: 


Ljv  =  mH  (—  -  1  )  g  (5) 

V  Pho  j 

This  shows  that  the  net  lift  is  independent  of  altitude.  At  the  pressure  altitude,  the  net  volume  becomes 
Vmax.  Let  the  density  ratio  at  this  altitude  be  denoted  by  ap.  The  equations  for  the  volume  at  the  pressure 
altitude  and  at  sea  level  are: 


Vmax  =  mH/  ( appH0 )  V0  =  mH/pH0  (6) 

It  follows  that  the  volume  of  Helium  at  sea  level  and  the  maximum  volume  are  related  through  ap. 

Vo  =  CTpVmax  (7) 

At  an  altitude  of  70  kft,  ap  =  0.061.  The  net  volume  at  sea  level  is  therefore  only  6%  of  the  total  available 
volume,  meaning  that  the  ballonets  must  occupy  the  remaining  94%.  It  will  likely  be  difficult  to  realize  this 
large  of  an  inflation  ratio.  This  problem  may  be  avoided  by  utilizing  vectored  thrust  and  aerodynamic  lift 
at  low  altitudes  for  climbing. 

To  maintain  vertical  equilibrium  through  the  buoyancy  force  alone,  the  net  lift  must  equal  the  combined 
weight  of  the  airship  structure,  systems  and  payload. 

Ln  =  (ms  +  mp)  g  (8) 

Here,  ms  represents  the  mass  of  the  structure  and  all  systems  (including  power,  propulsion,  pressure  control, 
etc.)  and  mp  represents  the  payload  mass. 


B.  Sizing  the  Airship  Hull 

Given  the  basic  aerostatic  principles  that  govern  the  airship,  and  assuming  a  classic  teardrop  shape  for  the 
hull,  it  is  possible  to  choose  a  set  of  dimensions  that  provide  a  volume  suitable  for  a  pressure  airship  operating 
at  70  kft.  We  first  combine  Eqs.  6,  5  and  8  to  find  the  following  expression: 


Vmax 


<-rp  ( Pao  Pho  ) 


(9) 


This  equation  determines  the  total  required  volume  of  the  airship  hull  based  upon  the  mass  of  the  airship 
structure  and  payload,  and  the  density  ratio  (which  is  a  function  of  altitude).  Because  the  mass  of  the 
airship  is  dependent  upon  the  total  volume,  it  becomes  a  design  challenge  to  maximize  the  volume  while 
minimizing  the  structural  mass.  Let  us  break  down  the  total  structural  mass  as  follows: 


ms  =  m.f  +  mo  (10) 

where  m/  is  the  mass  of  the  external  fabric  and  mo  is  the  remaining  mass  of  the  airship  structure  and  all 
supporting  systems.  The  mass  of  the  external  fabric  is  proportional  to  the  surface  area  through  the  fabric 
density,  pf. 

m.f  =  pfS  (11) 
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From  this  perspective,  a  spherical  shape  is  desired  as  it  provides  maximum  volume  with  minimum  surface 
area.  A  sphere  is  a  poor  choice  for  the  aerodynamics,  however,  as  it  generates  significant  pressure  drag.  The 
challenge  is  to  develop  a  shape  which  meets  the  competing  objectives  of  minimizing  both  the  drag  and  the 
surface  area. 

Recall  the  double  ellipsoid  geometry  of  the  hull  (see  Figure  1  on  page  3).  The  diameter  of  the  hull  is 
defined  as  D  =  26,  while  the  length  is  L  =  a\  +  a.2-  Let  us  now  introduce  n  as  the  length  ratio,  defined  as 
follows: 


(12) 


In  order  to  maintain  the  teardrop  shape,  the  constraint  of  k  >  1  is  imposed.  An  example  geometry  is  shown 
in  Figure  2,  with  L  =  250  m,  D  =  75  m,  and  n  =  2.  It  follows  that  ai  =  83.3  and  02  =  166.6. 


x [m]  y[m] 


Figure  2.  3D  View  of  Double-Ellipsoid  Hull  Geometry 


For  a  given  shape,  we  can  define  the  ratio  of  the  surface  area  over  the  volume  as  ksv.  The  equation  for 
maximum  required  volume  may  now  be  written  as: 


m0  +  Trip 


<?p  ( Pao  -  pH 0)  -  pfks 


(13) 


The  volume  of  this  geometry  is  736,311  m3,  and  the  surface  area  is  48,054  m2.  This  gives  a  value  of  0.0653 
for  ksv.  Using  this  geometry,  we  may  now  solve  for  mo  as  a  function  of  the  other  parameters: 


mo  —  Vmax  [dp  ( Pao  Pho )  Pfksv\  (14) 

This  represents  the  mass  available  for  the  internal  structure  and  all  systems.  Using  the  above  volume  and  a 
payload  mass  of  4,000  lb  (or  1,814  kg),  Figure  3  on  the  next  page  plots  several  curves  of  mo  as  it  varies  with 
the  fabric  density  and  the  surface-to-volume  ratio.  The  blue  dot  indicates  an  available  mass  of  19,651  kg 
corresponding  to  a  fabric  density  of  Pf  =  0.5  kg/m2  and  ksv  =  0.0653.  Polyester  has  traditionally  been 
the  preferred  fabric  for  airships,  and  has  a  density  of  0.33  kg/m2.  However,  additional  coatings  or  films  are 
required  to  reduce  permeability  and  increase  thermal  insulation,  which  adds  to  the  weight. 


C.  Drag  Estimation 

With  the  size  and  shape  of  the  airship  now  defined,  we  can  calculate  the  expected  drag  at  the  desired 
operating  condition.  The  usual  expression  for  the  aerodynamic  drag  force  on  a  body  is  given  as: 

D  =  1 pU2ACD  (15) 

where  p  is  the  atmospheric  density,  U  is  the  free-stream  velocity,  A  is  the  reference  area  and  Cd  is  the  non- 
dimensional  drag  coefficient  of  the  body.  The  drag  on  a  typical  airship  body  has  significant  contributions 
from  both  skin  friction  and  pressure. 
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Figure  3.  Mass  Available  for  Internal  Structure  and  Systems  vs.  Fabric  Density 


For  airships,  it  is  common  practice  to  express  the  reference  area  in  terms  of  the  hull  volume: 

A  =  V2/3  (16) 

The  drag  coefficient  associated  with  this  reference  area  definition  is  provided  in  2  for  several  values  of  thickness 
ratio.  The  thickness  ratio  is  the  maximum  hull  diameter  divided  by  the  total  length.  The  geometry  presented 
in  Section  B  on  page  4  has  a  reference  area  of  8,154  m2  and  a  thickness  ratio  of  0.3.  For  this  ratio,  the 
analytic  and  experimental  values  of  Cd  are  0.0244  and  0.0255,  respectively.  An  average  value  of  0.025  is 
chosen  for  subsequent  analysis.  If  we  instead  use  the  surface  area  S  as  the  reference  area,  the  corresponding 
drag  coefficient  is  0.00424. 

Consider  a  flight  condition  where  the  airship  is  at  its  target  altitude  of  70  kft,  and  flying  at  the  maximum 
expected  airspeed  of  46  m/s.  Using  a  drag  coefficient  of  0.025  and  the  given  reference  area,  the  total  drag 
force  is  approximately  15.4  kN.  Flying  at  the  nominal  airspeed  of  18  m/s,  the  drag  reduces  to  just  2.4  kN.The 
drag  will  be  higher  at  lower  altitudes,  but  we  are  interested  in  the  nominal  operating  condition,  where  the 
airship  spends  the  vast  majority  of  its  time. 

D.  Power  Requirements 

The  power  required  by  a  propulsion  system  with  efficiency  r]p  to  apply  thrust  T  at  velocity  U  is  given  as: 

TU 

Preq  =  —  (17) 

VP 

To  maintain  equilibrium,  the  thrust  must  equal  the  drag. 

Power  is  generated  by  the  solar  cells  during  the  day.  It  varies  from  zero  at  sunrise  to  a  maximum  value 
mid-day,  and  back  to  zero  at  sunset.  The  expression  is  given  as: 

Pgen{t)  =  SRIsoi(t)r)sc  (18) 

where  S  is  the  total  surface  area  of  the  hull,  R  is  the  percentage  of  the  surface  area  covered  with  solar  cells, 
Isoi  is  the  solar  irradiance,  and  rjsc  <  1  is  the  solar  cell  efficiency.  This  calculation  ignores  the  fact  that  the 
cell  normals  point  in  different  directions.  This  may  be  accounted  for  by  adjusting  r/ac  so  that  it  reflects  the 
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average  efficiency,  with  directional  effects  included.  The  total  energy  required  over  a  24  hour  period  is  found 
by  integrating  Preq  from  t  =  0  to  24  hours.  Using  an  average  velocity  U  and  substituting  the  drag  in  for  the 
thrust,  the  result  is: 

Ereq  =  +  Po)  tday  (19) 

where  Pq  represents  the  power  required  by  the  payload  (15  kW)  and  tday  is  the  number  of  seconds  in  a  day 
(left  as  a  parameter  for  later  manipulation).  The  total  energy  generated  over  a  24  hour  period  is  found  by 
integrating  Preq  from  sunrise  (tsr)  to  sunset  (fss).  The  variation  of  the  solar  irradiance  may  be  approximated 
as  a  sine  curve  from  0  to  180  deg.  The  result  is: 


Egen  —  S  RT)scEmax  (tss  tsr ) 

7T 


(20) 


where  Emax  is  the  peak  solar  irradiation  occurring  at  mid-day.  To  prevent  a  net  energy  loss,  the  amount 
generated  must  be  equal  to  or  greater  than  the  amount  required.  We  may  set  Eq.  19  equal  to  Eq.  20  and 
solve  for  the  minimum  required  solar  cell  coverage  ratio. 


7T  /  pUCD 

2  ^  2r]p 


-9l) 

SrjscJ 


1 


EmaxVscf s 


(21) 


where  fsun  is  the  fraction  of  the  day  receiving  sunlight,  a  parameter  that  varies  with  season  and  latitude. 

We  can  compute  the  required  solar  cell  coverage  ratio  R  for  a  range  of  efficiencies.  Let  fsun  =  0.5, 
Emax  =  500  W/m2,  and  Po  =  15  kW.  We  use  the  airship  geometry  from  Section  B  on  page  4,  and  choose  an 
altitude  of  70  kft.  Two  sets  of  curves  are  plotted  in  Figure  4,  each  set  corresponding  to  a  different  velocity. 
The  solid  curves  correspond  to  the  maximum  velocity  of  46  m/s,  while  the  dashed  curves  are  for  the  nominal 
velocity  of  18  m/s. 


propulsion  efficiency,  r| 


Figure  4.  Required  Solar  Cell  Coverage  Ratio  vs.  Propulsion  Efficiency 


Clearly,  R  is  physically  restricted  to  be  less  than  0.5,  since  the  solar  cells  must  be  located  on  the  top  half 
of  the  hull.  The  plots  show  that  a  feasible  coverage  ratio  is  possible  for  this  geometry  over  a  wide  range  of 
propulsion  and  solar  cell  efficiencies.  A  similar  analysis  is  presented  in  3. 

III.  Aerodynamic  Model 

Few  publications  are  available  which  address  the  aerodynamic  modeling  of  airships.  A  wealth  of  research 
was  performed  by  Maxwell  Munk  in  the  1920’s  and  1930’s,  which  provided  the  basis  for  subsequent  modeling 
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and  design  by  others.4,5  The  aerodynamic  model  presented  here  was  developed  using  the  procedure  outlined 
by  Jones  and  DeLaurier.6  This  model  includes  expressions  for  axial  force,  normal  force,  and  pitching  moment 
on  an  axisymmetric  airship  hull  with  4  equally  sized  tail  fins  -  2  horizontal  and  2  vertical.  The  equations 
are  valid  for  unseparated  flow  only. 

The  geometric  configuration  of  the  complete  airship  with  fins  and  gondola  is  shown  in  Figure  5.  It  is 
flying  at  an  air  speed  velocity  of  Vo  and  at  an  angle  of  attack  a.  The  sideslip  angle  (3  (not  shown)  is  defined 
as  a  positive  rotation  about  the  z-axis,  measured  from  the  x-axis.  The  forces  are  labeled  as  X,  Y,  Z ,  and 
the  moments  as  L,  M,  N.  The  forces  and  moments  on  the  hull  are  evaluated  from  the  nose  to  the  start  of 
the  fins.  Aft  of  this  point,  the  hull  and  fins  are  evaluated  together. 


Figure  5.  Steady-State  Aerodynamic  Model 


The  origin  of  the  body  frame  is  located  at  the  airship’s  center  of  volume.  For  the  double-ellipsoid 
geometry,  the  center  of  volume  is  located  on  the  x-axis  at  the  point: 

3 

xcv  =  ai  +  -  (a2  -  ai)  (22) 

O 

where  xcv  is  measured  backwards  from  the  nose.  For  the  airship  geometry  considered  in  this  analysis, 
xcv  =  114.583.  The  body  frame  must  be  located  at  the  CV  location  so  that  the  added  mass  and  inertia 
Let  p,  q,  and  r  denote  the  angular  velocities  about  the  x,  y,  and  z  axes,  respectively.  Let  Srjjdt  and 
Srudb  indicate  the  deflections  of  the  top  and  bottom  trailing  edge  flaps  of  the  rudder.  Similarly,  let  Selvl 
and  Selvr.  be  the  deflections  of  the  left  and  right  trailing  edge  flaps  of  the  elevator.  The  equations  for  the 
aerodynamic  forces  and  moments  are  provided  below: 

X  =  \pxi  cos2  a  cos2  f3  +  Cx2  sin(2a)  sin(a/2)]  (23) 

Y  =  i pV0 2  [CY\  cos(/?/2)  sin(2/3)  +  CY 2  sin(2 0)  +  CY 3  sin (/?)  sin(|/?|)  +  CY a{Srudt  +  SRUDb)\  (24) 

z  =  7;pVo  [Czi  cos(a/2)  sin(2a)  +  CZ2  sin(2a)  +  Cz 3  sin(a)  sin(|a|)  +  CZa(5Elvl  +  Selvr )]  (25) 

L  =  -pVg2  [Cli{Selvl  ~  Selvr  +  Srudb  -  Srudt)  +  Cl2  sin(/3)  sin(|/3|)]  (26) 

M  =  ipVp2  [CMi  cos(a/2)  sin(2a)  +  CM 2  sin(2a)  +  CM 3  sin(a)  sin(|a|)  +  CM4{Selvl  +  Selvr)]  (27) 

N  =  \pvo  icNi  cos(/3/2)  sin(2/3)  +  CN2  sin(2/3)  +  CN3  sin(/3)  sin(|/3|)  +  CN4(SElvl  +  Selvr)]  (28) 

(29) 

The  coefficients  are  direct  functions  of  the  airship  geometry.  They  are  fully  defined  in  the  appendix.  A 
complete  aerodynamic  model  is  available  in  the  Aircraft  Control  Toolboxa  for  Matlab. 

ahttp:  / /www.psatellite.com/products/html/ACT 
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IV.  Equations  of  Motion 


The  standard  equations  of  motion  for  a  rigid  body  with  6  degrees  of  freedom  are  used  to  model  the 
airship  dynamics.  In  addition,  the  effects  of  added  mass  and  inertia  are  included.  Beginning  with  the 
standard  6  DOF  equations,  we  have: 

mv  =  F  (30) 

=  T  (31) 

where  m  is  the  mass,  [/]  is  the  inertia  matrix  about  the  center  of  gravity  (CG),  v  is  the  linear  velocity  of 
the  CG,  uj  is  the  angular  velocity  of  the  body-fixed  frame,  F  is  the  sum  of  all  external  forces,  and  T  is  the 
sum  of  all  external  torques.  These  equations  are  valid  only  with  respect  to  an  inertial  reference  frame. 

It  is  convenient  to  express  the  linear  and  angular  velocity  of  aircraft  with  respect  to  a  body-fixed  frame. 
To  maintain  generality,  we  define  the  origin  of  the  body  frame  at  a  point  that  is  not  coincident  with  the 
CG.  This  is  useful  in  cases  where  the  CG  changes  slowly  during  flight.  It  also  enables  us  to  define  the  body 
frame  in  a  way  that  facilitates  the  inclusion  of  added  mass  and  inertia  terms  later.  Thus,  the  equations  of 
motion  are  expressed  in  a  non-inertial,  body-fixed  reference  frame  located  at  an  arbitrary  point  on  the  body. 
Figure  6  illustrates  the  inertial  frame  and  the  body-fixed  frame  on  the  airship.  Letting  [M]  denote  the 


Figure  6.  Coordinate  Frames 


transformation  matrix  from  the  inertial  to  body  frame,  the  position  of  the  CG  is  given  as: 

rc  =  r0  +  [M]TrG  (32) 

Assuming  that  derivative  of  rG  is  zero,  the  velocity  of  the  CG  is  expressed  in  rotating  body  frame  as  follows: 

Vc  =  V0  +  vG/0  (33) 

Vc  =  V0  +  uj  x  rG  (34) 

where  vG/0  indicates  the  velocity  of  the  CG  with  respect  to  the  origin  of  the  body  frame.  Differentiating 
this  velocity  with  respect  to  time,  the  translational  equation  of  motion  becomes: 

m  ( v'o  +  wXD0-rGXu  +  ux(wx  rG))  =  F0  (35) 


where  the  force  Fo  is  defined  in  the  body  frame.  The  corresponding  rotational  equation  of  motion  is: 


[/0]  uj  +  uj  x  ([I0]  uj)  +  mrG  x  (v0  +  uj  x  v0)  =  T 


(36) 


where  [To]  is  the  inertia  matrix  taken  about  the  origin  of  the  body  frame.  Combining  Eqs.  35  and  36,  we 
obtain  the  coupled  nonlinear  equations: 


rriE  —mr£, 
mr£  Iq 


Vo 

+ 

m  (uj  x  v'o  +  iv  x  (uj  x  rG )) 

'  Fo  ' 

C V 

uj  x  ([I0]  uj)  +  mrG  x  (uj  x  v0) 

T 

(37) 
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where  E  is  the  3x3  identity  matrix,  and  r£,  is  the  skew  symmetric  matrix  of  r^;,  defined  as: 


0 

—z 

y 

X 

z 

0 

— X 

with  r<3  = 

y 

.  ~ y 

X 

0 

z 

If  ra  =  0,  then  the  force  and  moment  equations  decouple  to  give: 


mE 

0 

V0 

+ 

m{uj  x  v0) 

F0 

0 

N 

UJ 

S' 

X 

c? 

£ 

T 

(38) 


(39) 


(40) 


which  are  the  conventional  6  DOF  equations  of  motion. 

Any  body  moving  through  a  fluid  displaces  that  fluid.  This  displacement  contributes  to  the  overall  change 
in  momentum  of  the  system.  When  considering  the  dynamics  of  conventional  aircraft,  the  momentum  change 
associated  with  the  displaced  air  may  be  ignored,  as  the  mass  of  displaced  air  is  extremely  small  compared 
to  the  aircraft.  For  an  airship,  however,  the  mass  of  the  air  that  it  displaces  is  roughly  equal  to  the  entire 
mass  of  the  vehicle.  Therefore,  it  is  necessary  to  account  for  the  linear  and  angular  momentum  of  the  air 
in  the  equations  of  motion.  This  results  in  “added  mass”  and  “added  inertia”  terms  added  to  the  nominal 
mass  and  inertia  of  the  airship. 

The  added  mass  and  inertia  associated  with  the  displaced  fluid  is  added  directly  to  the  nominal  mass 
and  inertia  on  the  left-hand  side  of  Eq.  37  on  the  page  before.  The  result  is  the  most  general  form  of  the 
6  DOF  equations  for  a  rigid  body  moving  through  a  fluid: 


mE  +  [_ M '] 

—mr£. 

Vo 

-1- 

m  (uj  x  v0  +  uj  x  (lu  x  rc)) 

F0 

mr£. 

io  +  ra 

uj 

uj  x  ([/0]  uj)  +  mrc  x  (uj  x  v0) 

T 

where  [. M ']  is  the  added  mass  matrix  and  [/,(]  is  the  added  inertia  matrix. 

The  sum  of  all  external  forces  and  torques  acting  on  the  airship  is  expressed  in  the  right-hand  side  of 
Eq.  41.  This  represents  the  aerodynamic  forces  and  moments,  the  buoyancy,  weight,  and  the  propulsive 
thrust  from  all  engines. 


A.  Computing  the  Added  Mass  and  Inertia  Matrices 


The  added  mass  and  inertia  matrices  are  functions  of  the  vehicle  shape.  For  the  double-ellipsoid  model  used 
in  this  analysis,  they  are  expressed  as  follows: 


ki 

0 

0 

’  0 

0 

0 

0 

0 

[Jo]  —  V(lair 

0 

CO 

0 

0 

0 

k2 

0 

0 

1 

CO 

(42) 


where  k\,  k2,  and  k2  are  ellipsoid  inertia  factors,  found  through  analysis  of  the  potential  flow  about  the 
ellipsoid.5  The  added  inertia  matrix  is  defined  with  respect  to  the  ellipsoid’s  center  of  volume  (see  Eq.  22  on 
page  8).  The  mass  of  displaced  air,  mair  ■  is  found  by  multiplying  the  volume  of  the  hull  by  the  local 
atmospheric  density,  where  the  volume  is  given  as: 

V  =  ^tt  («i  +  a2)  b2  (43) 

The  inertia  factors  for  a  single  axisymmetric  ellipsoid  are  pure  functions  of  the  semi-major  axis,  a, 
and  semi-minor  axis,  b.  A  complete  derivation  of  these  factors  may  be  found  in  a  1926  NACA  report  by 
Tuckerman.5  The  results  are  repeated  here  for  convenience.  First,  we  define  the  following  parameters: 


e  = 


(44) 
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/  =  log 


1  +  e 


1  —  e 

1  —  e2 


9  = 


a  =  2  g 


f 


—  e 


0  =  \-gi 

ez  2 

The  inertia  factors  for  this  mean  ellipsoid  are  then  expressed  as: 

a 


ki  =  - 

k2  =  - 


2  —  a 

P 

2-/3 
1 


k3  =  --  x 


( b 2  —  a2)  (a  —  (3) 


5  2  (62  —  a2)  +  (&2  +  a2)  (/3  —  a) 

Here,  we  define  a  as  the  mean  value  of  ai  and  a2  from  the  double-ellipsoid  hull  model. 


(45) 

(46) 

(47) 

(48) 

(49) 

(50) 

(51) 


V.  Open-Loop  Dynamics 

The  open-loop  dynamics  are  naturally  decoupled  into  two  modes:  longitudinal  and  lateral-directional. 
Separate  control  are  developed  for  each  mode.  The  longitudinal  mode  consists  of  the  following  states: 

•  q  -  Pitch  rate 

•  u  -  Axial  velocity 

•  w  -  Normal  velocity 

These  states  are  affected  by  the  following  three  controls: 

•  r  -  Throttle 

•  /r  -  Propeller  pitch  angle 

•  Selv  ~  Symmetric  elevator  deflection 

The  lateral-directional  mode  consists  of  the  following  states: 

•  p  -  Roll  rate 

•  r  -  Yaw  rate 

•  v  -  Lateral  velocity 

These  states  are  affected  by  only  two  controls: 

•  Sr ud  ^  Symmetric  rudder  deflection 

•  Sail  Differential  elevator  and  rudder  deflection  (termed  “aileron”  deflection  for  convenience) 

Linear  models  of  both  modes  are  obtained  directly  from  the  non-linear  model  by  application  of  the  central 
value  theorem.  The  velocity  components  of  the  linear  models  are  treated  as  perturbations  about  the  trim 
velocity.  Before  the  linearization  may  be  conducted,  the  trim  control  settings  must  be  computed.  Conditions 
for  trim  flight  are  discussed  in  Section  A  on  the  following  page. 
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Consider  the  airship  in  level  flight  at  21  km  altitude,  flying  at  an  airspeed  of  15  m/s.  The  linear  models 
at  this  flight  condition  are: 


’  0.0132 

-3.1  x  10-5 

-4.7  x  10"4  " 

’  0.0036 

1.1  x  104 

-0.0030  ’ 

■A-lon  — 

0 

-0.0042 

0 

Blon  — 

0.4737 

0.0139 

0 

7.079 

-5.9  x  10-4 

-0.0497 

0.0672 

-0.0208 

-0.0574 

(52) 


-1.6  x  10-9 

2.1  x  10"10 

-2.2  x  10-13  ’ 

0.0130 

0.0260  ’ 

A-lat  — 

-1.3  x  10-11 

0.0131 

4.7  x  10-4 

Blat  — 

0.0030 

0 

-3.0  x  10-9 

-7.079 

-0.0497 

-0.0574 

0 

This  model  corresponds  to  the  geometry  discussed  in  Section  B  on  page  4.  The  weight  is  set  equal  to  the 
buoyancy,  and  the  CG  is  coincident  with  the  CB. 

A.  Conditions  for  Trim  Flight 

The  airship  is  in  trim  flight  when  the  net  forces  and  moments  sum  to  zero.  This  includes  contributions  from 
the  aerodynamics,  gravity,  buoyancy,  and  thrust.  Three  control  effectors  are  used  to  establish  trim  flight: 
symmetric  elevator  deflection  ( Selv ),  throttle  (r  G  [0,1]),  and  propeller  pitch  angle  (//).  The  propeller 
pitch  angle  is  used  for  thrust-vectoring,  where  the  propellers  are  rotated  about  the  y-axis  of  the  airship.  We 
assume  two  identical  propellers  are  mounted  on  opposite  sides  of  the  gondola,  and  that  they  are  rotated 
together.  The  equations  for  longitudinal  trim  flight  are: 


EM* 

—  ZpTx  +  XpTy  +  qCM4$ELV  +  AIyD  +  MA^  —  o 

(53) 

EF* 

=  Tx  +  FAD  +  FAS  =  0 

(54) 

E* 

=  —  Tz  +  qCzA&ELV  +  FxD  +  Fxs 

(55) 

where  xp  and  zp  are  the  coordinates  of  each  propeller,  Tx  and  Tz  are  the  axial  and  normal  forces  from 
the  total  propeller  thrust,  and  q  is  the  dynamic  pressure.  The  aerodynamic  coefficients  Cm4  and  Cza  are 
defined  in  the  appendix.  The  superscripts  AD  and  AS  refer  to  aerodynamic  and  aerostatic,  respectively. 
The  aerodynamic  terms  are  computed  from  the  aerodynamic  model,  and  the  aerostatic  terms  are  found  to 
be: 


Mas 

y 

=  (B  -  mg)  [(zb  -  zg)  sin(0  +  a)  +  (xb  -  xg)  cos (9  +  a)] 

(56) 

F*S 

=  (B  —  mg)  sin(0  +  a) 

(57) 

Fas 

=  (B  —  mg)  cos(0  +  a) 

(58) 

where  B  is  the  buoyancy  force,  mg  is  the  weight,  <  Xb,Zb  >  are  the  coordinates  of  the  center  of  buoyancy 
(CB),  and  <  xq,zg  >  are  the  coordinates  of  the  center  of  gravity  (CG).  Clearly,  if  the  buoyancy  equals  the 
weight,  and  if  the  CG  and  CG  are  coincident,  then  there  exist  no  aerostatic  forces. 

Let  Tmax  be  the  maximum  achievable  thrust  for  each  propeller.  The  axial  and  normal  thrust  values  are 
then  defined  as: 

Tx  —  2r Tmax  cos  n  Tz  —  2TrIrnax  sin  //  (59) 

The  solution  for  trim  control  settings  is  readily  found  by  treating  Tx,  T,  and  Selv  as  the  unknowns  in 
Eqs.  56  to  58,  and  solving  the  3x3  matrix  inverse.  The  analytic  solution  is: 


Tx  = 

-(Fxad  +  Fas) 

(60) 

Tz  = 

kCzA  [zP  (FXD  +  Fas)  -  (MAD  +  Mas )]  +  kCM 4  ( FAD  +  Fas) 

(61) 

Selv  = 

yq  [*P  (F?D  +  Fxs)  -  (Fad  +  Fas)  (Mad  +  Mf)] 

(62) 

where  the  parameter  k  is 

defined  as: 

k  =  — - — — 

(63) 

Gm4  +  XpCz4 
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With  Tx  and  T,  known,  the  resulting  solutions  for  the  propeller  pitch  angle  and  throttle  are  found  from 
Eq.  59  on  the  preceding  page. 

Eq.  62  on  the  page  before  shows  that  the  magnitude  of  the  elevator  deflection  is  inversely  proportional  to 
the  dynamic  pressure,  q.  Since  the  aerodynamic  forces  and  moments  are  directly  proportional  to  q ,  it  cancels 
on  these  terms,  but  it  remains  on  the  aerostatic  terms.  As  the  velocity  is  decreased,  q  approaches  zero.  Thus, 
if  the  aerostatic  forces  and  moments  are  non-zero,  the  required  elevator  deflection  becomes  sensitive  to  the 
following  quantity: 


pAS  _  pAS  __  J^AS 

x  a,pi-  z  iviy 


(64) 


It  becomes  worthwhile  in  the  design  stage  to  minimize  this  quantity  through  careful  selection  of  the  CG,  CB 
and  propeller  location. 


VI.  Control  Law  Design 

An  integrated  guidance  and  control  system  has  been  designed  for  the  high  altitude  airship.  The  overall 
architecture  of  the  system  is  shown  in  Figure  7.  A  remote  operator  commands  the  target  latitude,  longitude 
and  altitude  for  the  airship.  This  information  is  provided  to  a  Trajectory  Planner  module,  which  also 
receives  the  measured  latitude,  longitude  and  altitude,  and  other  possible  data,  such  as  the  estimated  wind 
velocity.  The  Trajectory  Planner  periodically  performs  online  trajectory  optimization  to  generate  a  time- 
tagged  trajectory  sequence  for  the  control  law  to  follow.7  The  trajectory  consists  of  heading  angle  i/>,  flight 
path  angle  9,  and  absolute  velocity  V. 


Figure  7.  Control  System  Block  Diagram 


The  m  and  c  subscripts  distinguish  the  measured  and  commanded  signals. 

A  simple  control  law  design  provides  robust  feedback  control  of  the  airship’s  angular  rates  and  velocity. 
The  design  flight  condition  is  straight  and  level  flight  at  the  target  altitude  of  21  km,  and  at  the  nominal 
airspeed  of  18  m/s.  Based  upon  the  analysis  of  the  plant  input  matrix  for  each  mode,  separate  controllers 
were  designed  for  the  following  loops: 

•  Longitudinal  Mode: 

—  U  -tO-  T 
—  U  -to- 

—  q  -to-  Selv 

•  Lateral-Directional  Mode: 


-  r  -to-  6rud 
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—  V  -to-  5ruD 

-  p  -to-  Sail 


The  controller  for  each  loop  is  a  PID  with  rate  roll-off.  The  transfer  function  has  the  form: 


K(s) 


as2  +  bs  +  c 
s(s  +  0.1) 


(65) 


where  a,  b ,  and  c  are  directly  related  to  the  proportional,  integral  and  derivative  terms  of  a  standard  PID 
controller.  These  coefficients  were  chosen  for  each  loop  through  an  iterative  process,  in  which  the  robustness 
margins  and  the  time  response  were  evaluated,  until  the  desired  performance  was  achieved.  The  standard 
military  robustness  requirements  (45  deg  phase  margin,  6  dB  gain  margin)  are  met  in  each  loop. 

A  simulation  of  the  complete  non-linear  model  is  conducted  with  the  controllers  implemented  at  1  Hz. 
Perfect  measurements  and  actuation  are  assumed.  The  airship  is  initialized  at  the  design  altitude,  flying 
at  24  m/s  with  a  2  degree  angle-of-attack  and  a  1  -degree  sideslip.  The  desired  velocity  is  25  m/s,  with 
zero  angle-of-attack  and  sideslip.  The  time-response  is  shown  in  Figure  8.  Even  though  this  airspeed  is  far 
from  the  speed  for  which  the  controllers  were  designed,  they  still  provide  good  performance.  The  augmented 
response  is  much  faster  than  that  of  the  open-loop  system,  and  all  of  the  actuators  stay  well  within  their 
physical  limits. 


Time  [sec] 


Figure  8.  Time  Response  of  Closed-Loop  System 


It  is  interesting  to  examine  the  closed-loop  stability  afforded  by  this  single  control  law  across  a  range 
of  flight  conditions.  In  this  case,  we  vary  the  airspeed  from  1  to  25  m/s,  keeping  the  altitude  constant  at 
21.3  km.  The  root  locus  for  the  both  the  longitudinal  and  lateral-directional  modes  is  shown  in  Figure  9  on 
the  following  page.  The  direction  of  increasing  velocity  is  denoted  by  increasing  marker  size.  In  each  system, 
the  dominant  closed-loop  poles  approach  the  imaginary  axis  as  the  velocity  decreases,  signifying  a  reduction 
in  stability. 

The  airship  design  considered  in  this  analysis  is  statically  neutral.  The  buoyancy  equals  the  weight,  and 
the  two  forces  act  at  the  same  point.  So,  in  the  absence  of  aerodynamic  forces,  the  vehicle  will  simply  hover. 
This  agrees  with  the  pole  locations  approaching  the  imaginary  axis  at  zero  velocity.  If  the  center  of  buoyancy 
is  different  from  the  center  of  gravity,  which  is  likely  in  a  real  system,  the  vehicle  would  be  statically  unstable. 
A  minimum  amount  of  airspeed  would  be  required  in  order  to  maintain  stability. 
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Longitudinal  System 
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Lateral-Directional  System 


Figure  9.  Root  Locus  of  Closed-Loop  Systems  at  Various  Airspeeds 


VII.  Conclusions 

This  paper  has  provided  useful  methods  for  the  modeling  and  control  of  a  typical  airship  configuration, 
and  discussed  the  unique  issues  associated  with  the  operation  of  airships  at  high  altitudes.  A  series  of  steps 
for  sizing  the  airship  hull  were  presented,  based  on  the  principles  of  aerostatics  and  daily  energy  balance 
constraints.  A  steady-state  aerodynamic  model  was  provided,  based  on  the  model  by  Jones.6  The  dynamics 
of  the  airship  in  flight  was  then  expressed  using  the  six  degree  of  freedom  equations  of  a  rigid  body  with  an 
offset  CG,  and  including  the  effects  of  added  mass  and  inertia.  Linearized  models  of  the  open-loop  dynamics 
were  then  used  to  develop  classical,  SISO  controllers  for  tracking  desired  angular  rates  and  linear  velocity. 
The  approach  has  shown  that  useful  models  of  airship  dynamics  may  be  readily  developed  and  analyzed. 
For  example,  the  model  developed  in  this  paper  has  been  used  in  the  development  of  trajectory  optimization 
algorithms. ' 

Future  research  on  this  topic  should  address  the  issue  of  model  validation.  High-fidelity  models  of 
airships  may  be  obtained  with  the  aid  of  CFD  tools  or  wind-tunnel  testing.  However,  these  methods  are 
prone  to  high  cost  and  long  development  cycles.  If  the  dynamics  can  be  modeled  with  sufficient  accuracy 
using  an  analytic  approach,  then  multiple  designs  and  analyses  may  be  carried  out  rapidly  at  low  cost.  A 
particular  issue  important  to  high-altitude  airships  is  the  extensive  flight  envelope  that  they  must  traverse 
altitudes  range  from  sea-level  to  over  20  km,  and  airspeed  velocities  fluctuate  strongly  with  the  wind.  It 
is  therefore  important  to  understand  the  dynamics  across  this  entire  envelope,  and  to  assess  the  robustness 
and  performance  capabilities  with  respect  to  wind.  It  appears  that  this  problem  lends  itself  well  to  modern, 
robust  control  methods.  The  development  of  MIMO  controllers  using  and  /x-synthesis  techniques,  for 
example,  would  offer  robustness  guarantees  in  the  face  of  uncertain  dynamics. 

Appendix 

The  force  and  moment  equations  defined  in  Section  III  on  page  7  contain  20  aerodynamic  coefficients. 
Each  of  these  coefficients  is  defined  below: 


Cxi 

Cx2 

Gy  2 

Gy3 

Cy4 


-  [ CDh0Sh  +  Cnf0Sf  +  CDg0Sg] 

(k2  -  k\ )  rjkhSh 


Gy  i  =  Czi 

_i  (dcL 


2  \  da  J  f 
G Z2 


sfVf 


-  [ CDchJlSh  +  CocfSf  +  CDcgSg \ 


(66) 

(67) 

(68) 

(69) 

(70) 

(71) 

(72) 
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(73) 

(74) 


Cz3 


Cz4 

—  [CDchJlSh  +  CocfSf ] 


C,A 

CL2 

Cmi 

Cm2 

Cm3 
Cm  4 

CNj 


C Deg  Sg  Igz 

—  (ki  -  k2)r)khShL 

—  [CDch,J2ShL  +  CocfSflf^l 

—  CMj 


(75) 

(76) 

(77) 

(78) 

(79) 

(80) 
(81) 


The  following  table  defines  all  of  the  parameters  appearing  in  the  coefficient  equations  (Eq.  66  on  the 
preceding  page  to  Eq.  81). 


Table  2.  Definition  of  Parameters  in  the  Aerodynamic  Coefficient  Equations 


Parameter 

Definition 

Value 

CDh0 

Hull  zero-incidence  drag  coefficient 

0.025 

CDf0 

Fin  zero- incidence  drag  coefficient 

0.006 

cDgo 

Gondola  zero-incidence  drag  coefficient 

0.01 

CDch 

Hull  cross-flow  drag  coefficient 

0.5 

CDcf 

Fins  cross-flow  drag  coefficient 

1.0 

C Dcq 

Gondola  cross-flow  drag  coefficient 

1.0 

(OCl\ 

\  dot  Jf 

Derivative  of  Hn  lift-coefficient  with  respect  to 

the  angle-of-attack  at  zero  incidence 

5.73 

(QCl\ 
l  95  If 

Derivative  of  Hn  lift-coefficient  with  respect  to 

the  flap  deflection  angle 

1.24 

sh 

Hull  reference  area,  F2/3 

8154  m2 

Sf 

Fin  reference  area 

3656  m2 

Sa 

Gondola  reference  area 

202  m2 

Ifl 

x-distance  from  origin  to  aerodynamic  center  of  fins 

117.5  m 

lf2 

x-distance  from  origin  to  geometric  center  of  fins 

129.7  m 

l.f3 

y ,  ^-distance  from  origin  to  aerodynamic  center  of  fins 

18.3  m 

Igx 

x-distance  from  origin  to  aerodynamic  center  of  fins 

29.2  m 

Igz 

^-distance  from  origin  to  aerodynamic  center  of  fins 

40.0  m 

V.f 

Fin  efficiency  factor  accounting  for  the  effect  of  the  hull  on  the  fins 

0.29 

Vk 

Hull  efficiency  factor  accounting  for  the  effect  of  the  fins  on  the  hull 

1.19 

h 

0.33 

h 

1_  ra2  xdA(x) 
ohih  J —a  1  ax 

-0.69 

Jl 

£  /-I  2r(x)dx 

1.31 

J2 

snz  f—n.i  2 r(x)xdx 

0.53 
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The  lift  and  drag  coefficients  were  approximated  based  upon  the  shape  of  the  various  components.  Data 
from  the  NACA  1412  airfoil  was  used  for  the  fins.  The  gondola  is  assumed  to  be  stream-lined,  with  a 
cross-sectional  area  larger  than  its  frontal  area.  See  Figure  5  on  page  8  for  an  illustration  of  the  geometry. 


Using  the  double  ellipsoid  geometry,  the  solution  to  the  integrals  are: 

=  '^(W2)  <82> 

h  =  (ai  -  2 a2f  -  3 a,/2}  -  tj-h  (83) 

Ji  =  3(^73  1  \  '  f-  f  I  2a2sin~1(/))  (84) 

J2  =  ■'1£4f=  +  3Z^7s(»2-«;-«l(l-/2)3/2)  (85) 


where  f  =  (lh  —  «i)/a2>  hi  is  the  length  of  the  hull  up  to  the  leading  edge  of  the  fins,  and  xcv  is  defined  in 
Eq.  22  on  page  8. 
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